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ABSTRACT 

The primary impediment to continued improvement of traditional charge-based electronic devices in accordance 
with Moore’s law is the excessive energy dissipation that takes place in the devices during switching of bits. One 
very promising solution is to utilize strain-mediated multiferroic composites, i.e., a magnetostrictive nanomagnet 
strain-coupled to a piezoelectric layer, where the magnetization can be switched between its two stable states 
in sub-nanosecond delay while expending a minuscule amount of energy of attojoule at room-temperature. 
Apart from devising digital memory and logic, these multiferroic devices can be also utilized for analog signal 
processing, e.g., voltage amplifier. First, we briefly review the recent advances on multiferroic straintronic devices 
and then we show here that in a magnetostrictive nanomagnet, it is possible to achieve the so-called Landauer 
limit (or the ultimate limit) of energy dissipation of amount kT ln(2) compensating the entropy loss, thereby 
linking information and thermodynamics. 

Keywords: Nanoelectronics, spintronics, multiferroics, magnetoelectrics, straintronics, analog signal processing, 
Landauer limit, information, thermodynamics, energy-efficient computing 

1. INTRODUCTION 

Electric field induced magnetization switching in multiferroic composites works according to the principles of 
converse magnetoelectric effect, i.e., when an electric field is applied across the structure, it induces a magnetic 
anisotropy, which can rotate the magnetization. This eliminates the need to utilize the cumbersome magnetic field 
or high spin-polarized current to rotate the magnetization.®^ Hence, it can harness an energy-efficient binary 
switch replacing the traditional charge-based transistors for our future information processing systems.®® This 
turns out to be a very promising mechanism since a small voltage can generate a huge magnetic anisotropy, e.g., 
with suitable choice of materials, when a voltage of few millivolts is applied across a strain-mediated multiferroic 
composite device (see Fig. [T]), i.e., a magnetostrictive nanomagnet attached to a piezoelectric layer,®® the 
piezoelectric layer is strained and the generated strain is transferred to the magnetostrictive layer. Then the 
induced stress anisotropy in the nanomagnet can switch the magnetization between its two stable states that 
store a binary information 0 or 1.13® This study has opened up a new field named straintronicJ^EnEII and 
experimental efforts to realize such devices are considerably emerging .E3III1 Although the experimental efforts 
have demonstrated the induced stress anisotropy in the magnetostrictive nanomagnets, the direct experimental 
demonstration of switching delay (rather than ferromagnetic resonance experiments^ and utilizing low-thickness 
piezoelectric layers [e.g., < 100 nm of lead magnesium niobate-lead titanate (PMN-PT)P^ for ultra-low energy 
dissipation are still under investigation. 

While utilizing strain-mediated multiferroic composites for digital binary switching of the magnetization 
has been investigated, the analog signal processing capability using inherently digital nanomagnets has not 
been thought of. It is true that computing and signal processing tasks are mostly shifted to digital domain, 
however, sometimes analog signal processing is fundamentally necessary, e.g., processing of natural signals.173 
When a transmitted signal is received at the receiver end, the signal is weak in magnitude and also noisy due 
to attenuation and noise in the environment. Hence, the signal needs to be amplified and filtered. Only after 
converting the analog signal to digital domain, we can process the signal digitally. Therefore, for different 
purposes, we do have the requirement of analog signal processing. In a very recent development, it is shown 
that we can utilize these multiferroic straintronic devices for ultra-low-energy analog signal processing too.^ A 
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Figure 1. Schematic diagram of a strain-mediated multiferroic composite (piezoelectric-magnetostrictive 
heterostructure), and axis assignment. The magnetostrictive nanomagnet is shaped like an elliptical cylinder and 
it has a single ferromagnetic domain. It’s mutually anti-parallel magnetization states along the z-axis store the binary 
information 0 or 1. In standard spherical coordinate system, 9 is the polar angle and <j) is the azimuthal angle. According 
to the convention, we term the a-axis the easy axis, the y-axis the in-plane hard axis, and the r-axis the out-of-plane hard 
axis based on the chosen dimensions of the nanomagnet. 


voltage applied across a multiferroic composite induces stress on the magnetostrictive layer and this modulates 
the potential landscape of the nanomagnet. To be able to harness the analog nature, we need to rotate the 
magnetization continuously rather than having an abrupt switching as required by digital computation. The 
continuous rotation of magnetization is conceived by considering that magnetization is not exactly stable at a 
position rather it is fluctuating around a mean value due to thermal agitations. The applied voltage modulates the 
potential landscape of the magnetostrictive nanomagnet and the mean value of magnetization can be changed 
continuously. Using tunneling magnetoresistance (TMR) measurement this continuous mean value can 

be compared with respect to a fixed nanomagnet and a continuous output voltage while varying the input 
voltage can be produced. Stochastic Landau-Lifshitz-Gilbert (LLG) equation of magnetization dynamics at 
room-temperatur^SMMI jg solved to demonstrate this concept of analog signal processing functionality.!^ 

1.1 Theory and model 

It is a common practice to model nanomagnets as shape-anisotropic elliptical cylinders (see Fig.[T]) to produce an 
energy barrier between its two stable states anti-parallel to each other along the elongated axis (x-axis). When 
stress is generated on the magnetostrictive nanomagnet, the total potential energy of the stressed polycrystalline 
(we neglect the magnetocrystalline anisotropy) single-domain nanomagnet per unit volume is the sum of the 
shape anisotropy energy and the stress anisotropy energy!^! 

Etotal{9, <(>, t) = EghapeiO, (j)) E E stress i9,t) 

— Eghapei^^^ stn 9 cos 9^ (1) 

where 

BshapeW = (1/2) M [Hk + Hd cos^cj}], (2) 

Bstressit) = (3/2) AsCr(t), (3) 

M = imqMs, fiQ is permeability of free space. Mg is the saturation magnetization, Hk = {Ndy — Ndz) Mg is 
the Stoner-Wohlfarth switching field,!^ Hd = {Ndx — Ndy) Mg is the out-of-plane demagnetization field,!^ Ndm 
is the m*^ (m = x,y,z) component of the demagnetization factor!^ {Ndx ^ Ndy > Ndz), (3/2)As is the 
magnetostriction coefficient of the magnetostrictive nanomagnet ct = Ue is the stress on the nanomagnet, Y 
is the Young’s modulus, e is the strain that generates the stress. 

Materials having positive As (e.g., Terfenol-D, Galfenol) requires a compressive (negative by convention) 
stress to favor the alignment of magnetization along the minor axis (y-axis in Fig.[T]), while the materials having 
negative As (e.g., Fe, Co, Ni) requires a tensile (positive by convention) stress for the same. The gradient 
of potential profile due to shape and stress anisotropy exerts an effective held on the magnetization and this 
generates a torque Te,!^ while the torque due to thermal huctuations Tth is treated via a random magnetic 
field.El^ 










In the macrospin approximation, the magnetization M of the nanomagnet is constant in magnitude but it 
varies with direction, so that we can represent it by a vector of unit norm rim = M/|M|. The magnetization 
dynamics under the action of these two torques Te and Tth is described by the stochastic Landau-Lifshitz- 
Gilbert (LLG) equatioiP^^f^ as follows: 

drim I dnm\ I 7 I r , 

where a is the phenomenological Gilbert damping parameter!^ through which magnetization relaxes to the 
minimum energy position, and 7 is the gyromagnetic ratio of electrons. Note that M, Te, and Txh are all pro¬ 
portional to the nanomagnet’s volume. Solving the above equation, we can track the trajectory of magnetization 
over time. 

For the magnetostrictive layer, we need to choose a material that maximizes the product (3/2)As Y. Terfenol- 
D (TbDyFe), which has 30 times higher magnetostriction coefficient in magnitude than the common ferromagnetic 
materials (e.g., Fe, Go, Ni), has the highest (3/2)Asy.l^ If it needs to avoid the rare-earth materials (e.g., Tb 
and Dy in Terfenol-D), we can also utilize Galfenol (FeGa),l^l^ which has 6 times less (3/2)As, but twice high 
Y than that of Terfenol-D. 1^1 

For the piezoelectric layer, we may use lead-zirconate-titanate (PZT), but using lead magnesium niobate-lead 
titanate (PMN-PT) is preferable since it has high piezoelectric coefficient and it can generate anisotropic strain, 
which allows us to work with lower voltage for a required strain, thereby reducing the energy dissipation. PMN- 
PT layer has a dielectric constant of 1000, (i3i=-3000 pm/V, and (i32=1000 pm/V.^^^With the piezoelectric layer’s 
thickness tpiezo=24 nm,l^ V = 1.9 mVs (2.9 mVs) of voltages would generate 20 MPa (30 MPa) compressive 
stress [a = Y def fiV/tpiezo), where deff = (dsi — d32)/(l + v)] in the magnetostrictive Terfenol-D layer, which 
has y = 80 GPaP and Poisson’s ratio v = 0.3.®l Modeling the piezoelectric layer as a parallel plate capacitor 
(~100 nm lateral dimensions), the capacitance C=2.6 fF and thus CV^ energy dissipation turns out to be < O.I 
aJ. This is the basis of ultra-low-energy computing using these multiferroic devices 

We do consider the distribution (rather than a fixed value) of initial orientation of magnetization due to 
thermal fluctuations, which has crucial consequence on device performance.^^^ Note that we have assumed 
here uniform magnetization and uniform strain in the magnetostrictive nanomagnet. This is quite valid in 
small length scales and such assumptions allow us to get critical insights on the device operation. A detailed 
space-dependent solution may be sought for quantitative purposes, however, it would require an enormous time 
to execute. Hence, we stick to the aforesaid assumptions. 

1.2 Digital memory 

Here we briefly review the binary switching of magnetization in the magnetostrictive nanomagnets for devising 
digital memory devices.^ The usual perception is that stress can rotate magnetization only by 90° from ±z-axis 
to Ty-axis (see Fig. HI). However, the torque due to stress acts along the out of magnet’s plane (e^ direction) and 
therefore magnetization lifts out-of-plane (although the out-of-plane demagnetization field due to small thickness 
of nanomagnet tries to keep the magnetization in-plane). 

Te, stress ~ X V Egi’pQss = (3/2) \g(jsiTi{f20^ (5) 

This out-of-plane excursion of magnetization generates an intrinsic asymmetry, which can completely switch 
the magnetization by 180°.!^ Full 180° switching is desirable since it facilitates having the full tunneling magne¬ 
toresistance (TMR) while electrically reading the magnetization state using a magnetic tunnel junction (MTJ). 
As depicted in the Fig. HJ magnetization’s initial distribution due to thermal fluctuations is quite wide.l^ Also, 
thermal fluctuations create a distribution for the time required to reach at the hard plane {9 = 90°, y-z plane). 
Therefore, to tackle this wide distribution, it requires a sensing circuitry to detect when magnetization reaches 
around 9 = 90° and withdraw/reverse stress subsequently to be able to switch 180° successfully.^^ The simu¬ 
lation results show that fast ramp rates and high stresses are conducive to successful switching.l^IHl jt should 
be noted that such switching leads to a toggle memory unless we have a mechanism to maintain the direction 
of switching .1^1^ According to Ref. [311 the interface coupling between the polarization and magnetization can 
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Figure 2. Digital memory using multiferroic straintronic devices. Switching of magnetization between two anti¬ 
parallel states (iz-axis). Stress rotates magnetization out of magnet’s plane and when magnetization reaches the hard 
plane {9 = 90°), intrinsic dynamics dictates that magnetization lifts out-of-plane in a certain direction so that a complete 
180° switching of magnetization is possible.® While switching along the -^-axis rather than -fy-axis (shown by arrows), 
the directions of the out-of-plane excursions would be exactly opposite.® Three distributions are shown: the one at 
6 — 180° depicts the initial distribution when no stress is active, the other two distributions around 9 — 90° and 6^ = 0° 
correspond to 60 ps ramp period and 15 MPa stress.® 

maintain the direction of switching. Also the strong coupling facilitates error-resilient switching without the 
need of having a sensing circuitry and it can lead to lowering the lateral dimensions of the nanomagnet to ^10 
nm. The 90° switching mechanisms can be also utilized to direct switching in a particular direction, however, it 
gives us lower TMR. 

On reading the magnetization state of the magnetostrictive nanomagnet, a material issue crops up since 
magnetostrictive materials (Terfenol-D, Galfenol) cannot be in general used as the free layer of an MTJ. Usually, 
CoFeB is used as the free layer alongwith MgO spacer that gives us high TMR,I^ while using half-metals can lead 
to even higher TMR.I^ This issue can be simply solved by magnetically couple the magnetostrictive layer and 
the free layer, e.g., exchange coupling or introducing an insulator between the layers and exploiting the dipole 
coupling in between to rotate them concomitantly.^^ 

1.3 Digital logic 

Apart from devising memory bits, these straintronic multiferroic devices have been proposed for logic design 
purposes too.^^^^ Figure [3^ depicts Bennett clocking mechanisirPS of information propagation through a 
horizontal chain of multiferroic straintronic devices.^ There is no wired connection between the nanomagnets 
and the information transfer happens due to dipole coupling between the nanomagnets. While considering room- 
temperature operation, the intrinsic dynamics due to out-of-plane excursion of magnetization causes switching 
failureJ^and it can remedied by using the intrinsic asymmetry, i.e., releasing/reversing stress dynamically when 
magnetization reaches the hard plane {0 = 180°)fS^or exploiting the asymmetry due to interface coupling.l^ 

An intriguing mechanism to build logic is to use the concept of read-write units.^ The written bit in the 
magnetostrictive nanomagnet is read by an MTJ (read-unit) and the output can be fed to next write units 
(multiferroic composites), i.e., individual read-write units can be concatenated since they have voltage gain 
and input-output isolation.!!^ Also, by utilizing two inputs as shown in the Fig. [SId, 2-input universal logic gates 
(NAND and NOR) can be designed.^ Higher input logic gates and majority logic gates can be devised according 
to the same concept.^ This design methodology can overwhelmingly simplify the design of a large scale circuit 
and portend a highly dense yet an ultra-low-energy computing paradigm. 

1.4 Analog signal processing 

While a bistable double-well separated by an energy barrier is desired for digital computing, analog signal 
processing requires a continuous rotation of magnetization. The potential landscape of a shape-anisotropic 
single-domain nanomagnet is a symmetric double-well and hence it needs to stop the abrupt switching between 
the two stable states for analog signal processing. We can make the potential landscape monostable by utilizing 
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Figure 3. Digital logic design using multiferroic straintronic devices, (a) Information propagation through a 
horizontal chain of straintronic devices uni-directionally. Since the dipole coupling between the neighboring nanomagnets 
is bidirectional, we need to impose the unidirectionality in time (using a 3-phase clocking scheme) to subsequently clock 
the nanomagnets in the chain.^^ With this scheme, logic gates like NAND and NOR, and complex circuits can be devised. 
(Reprinted with permission from Ref. 1111 Copyright 2014, AIP Publishing LLC.) (b) Logic design using read-write 
units. The read unit is an MTJ, which reads the free layer’s magnetization direction with respect to the hxed layer. 
The magnetostrictive layer and the free layer are dipole coupled (anti-parallel orientation), hence the magnetostrictive 
nanomagnet’s direction can be read. Such read-write units can be concatenated to devise complex circuits. Two inputs 
can be incorporated to design single-element straintronic universal logic gates (e.g., NAND, NOR).^By applying voltages 
at the terminals A and B, the magnetization of the magnetostrictive nanomagnet can be switched and depending on the 
design of the nanomagnets, the respective logic operation can be performed.!^ 
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Figure 4. Continuous rotation of magnetization in magnetostrictive nanomagnets considering thermal fluc- 
tnations. (a) Magnetization is fluctuating due to room-temperature (300 K) thermal agitations around one easy axis 
9 — 180°. The distribution is achieved by solving stochastic LLG at 300 K. The distribution turns out to be a Boltzmann 
distribution as expected, (b) The potential landscape of the magnetostrictive nanomagnet is modulated with stress. As 
the energy barrier decreases, the mean orientation of the magnetization changes continuously rather than abruptly. When 
sufficient stress is exerted, magnetization comes toward the hard-axis {6 — 90°). Note that the potential landscape of 
the magnetostrictive nanomagnet is made monostable (while no stress is active) by having a magnetic coupling from a 
neighboring nanomagnet to avoid abrupt switching. 
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Figure 5. Design of a voltage amplifier using multiferroic straintronic devices, (a) The device structure to harness 
the analog nature consists of a multiferroic composite attached to an MTJ separated by an insulator. The magnetostrictive 
nanomagnet is magnetically coupled to the free layer via dipole coupling. Iread is a constant read current and the output 
voltage Vout varies continuously with Vin- (b) The relation between the input voltage Vin and the output voltage Vout for 
a Terfenol-D/PMN-PT multiferroic composite. This is an inverter like characteristics with a high-gain transition region. 
If we bias the input voltage in this high-gain region, we would get an amplified output voltage. 


magnetic coupling from a neighboring nanomagnet. As depicted in Fig. 21 we choose the potential well at 
9 = 180°. Figure UK shows magnetization’s fluctuations due to room-temperature (300 K) thermal agitation 
by solving stochastic LLG and it fits quite well with the Boltzmann distribution as expected. Figure Hb shows 
that with the application of stress, magnetization’s mean orientation shifts continuously and when the stress 
is sufficient enough, magnetization comes close to 9 = 90°. An amplifier with a voltage gain can be designed 
utilizing this methodology as depicted in the Fig. [5l 

We can directly put an AC signal at the input terminal of the device (see Fig. [5b), solve stochastic LLG, and 
the output voltage can be extracted from the TMR measurement of the MTJ. For an AC signal with 1 GHz 
frequency, a voltage gain of 50 is achieved, while dissipating a miniscule amount of energy of ~0.1 attojoule/cycle 
at room-temperature.^ A rectifier characteristics can be obtained from these straintronic devices too since with 
the application of a voltage of opposite polarity (that generates stress of opposite polarity too, see Fig. |4b), 
magnetization becomes more confined in the 9 = 180° well. The degree of rectification capability increases with 
the TMR of the MTJ.ESl 


2. LANDAUER LIMIT OF ENERGY DISSIPATION 


Computers are physical systems and therefore a computation or processing of information is subjected to physical 
principle, e.g., the second law of thermodynamics. According to Landauer’s principle,^ in a classical system with 
two degenerate ground states, a minimum amount of energy proportional to temperature kTln(2) (~ 3 x 10“^^ 
joule at T=300 K, k is the Boltzmann constant, and T is temperature) must be dissipated to erase a bit of 
information compensating the entropy loss, thereby linking the information and thermodynamics. Maxwell 
introduced his controversial demor^^^^ over a century ago elucidating the relationship between information and 
entropy, however, the famous demon was later exorcized by Bennett Note that it is possible to have 
stochastic violation of Landauer’s principle in small systems, since they are prone to thermal fluctuations.^^^^^ 
Stochastic violations of second law have been experimentally observed,^^ however, the second law is still 
safeguarded on averageW^ 


Recently, Landauer limit of energy dissipation is experimentally demonstrated for a colloidal particle with 
linear motion.l^ However, the Landauer’s limit is not experimented for a rotational body like the magnetiza¬ 
tion in a nanomagnetWith the advent of experimental apparatus, single-domain nanomagnet^^IlLsHllol 












Figure 6. Deterministic switching of magnetization to a final stable state in a magnetostrictive nanomagnet. 

(a) A shape-anisotropic single-domain magnetostrictive nanomagnet having two stable states along the elongated direction 
and axis assignment, (b) The time-cycle of uniaxial stress (along 2 -direction) and an asymmetry making field, H (along 
0 = 0° direction) acting on the magnetostrictive nanomagnet, (c) The complete three-dimensional potential landscapes 
of the nanomagnet at different instants of time. The solid lines on the landscapes correspond to the in-plane {y-z plane, 
<j) = ±90°) potential profiles. Although the in-plane potential landscapes are the minimum energy positions, the out-of- 
plane potential landscapes do play an important role in shaping the magnetization dynamics as discussed in the text. At 
time ti, magnetization can reside in either of the wells (0 = 0° or 0 = 180°) having 50% probability each, however, at time 
t 2 , it reaches at 0 = 90°, i.e., the bit is erased. At time ts and onwards, magnetization reaches at 0 ~ 0° deterministically 
due to the asymmetry making field. 


having two stable states separated by an energy barrier can become the staple of future information processing 
systems (since the excessive energy dissipation during switching of bits has been the bottleneck behind utilizing 
traditional charge-based transistor electronics furtheJ^i^^^ . Particularly, magnetostrictive nanomagnets in mul- 
tiferroic heterostructures have profound potential to act as the basic building block in ultra-low-energy computing 
systems.l^li^ Here, we show that Landauer limit of energy dissipation is achievable in magnetostrictive nanomag¬ 
nets. It is pointed out that magnetization may deflect out of magnet’s plane during its dynamical motion and 
even a very small out-of-plane excursion plays a crucial role in shaping the magnetization dynamics. Therefore, it 
is imperative to consider the complete three-dimensional potential landscape and full three-dimensional motion 
of the magnetization, rather than assuming an overdamped particle with linear motion.1^1^ We have solved 
stochastic Landau-Lifshitz-Gilbert (LLG) equatiorP^HH! jn the presence of thermal fluctuations to determine the 
energy dissipation during the erasure of a bit of information and show that stochastic violation of the Landauer 
bound is possible, nonetheless the bound is respected on average. 

2.1 Model 

Figure|6^ shows a single-domain magnetostrictive nanomagnet shaped like an elliptical cylinder with its elliptical 
cross-section lying on the y-z plane; the major axis and minor axis are aligned along the 2 - and y-direction, 
respectively. In standard spherical coordinate system, 6 is the polar angle and (j) is the azimuthal angle. The 
magnet’s plane is y-z plane (</> = ±90°). Any deviation of magnetization from magnet’s plane is termed as 
out-of-plane excursion. The dimensions of the major axis, the minor axis, and the thickness are a, 6, and I, 
respectively. So the magnet’s volume is H = (7r/4)a&L Due to shape anisotropy, the two anti-parallel degenerate 
states 0 = 0° and 6 = 180° along the 2 -direction (easy axis) can store a binary bit of information 0 or 1. The 
y-axis and x-axis are the in-plane and the out-of-plane hard axis, respectively. Since I b, the out-of-plane hard 
axis is much harder than the in-plane hard axis. 








Figureinb shows the time-cycle of uniaxial stress (along z-direction) and an asymmetry making field H (along 
0 = 0° direction) acting on the magnetostrictive nanomagnet. The adiabatic pulses should be slow enough that 
magnetization follows its potential landscape and the quasistatic assumption is valid. We can write the total 
energy of the nanomagnet per unit volume as the sum of three energies - the shape anisotropy energy, anisotropies 
induced due to applied stress and asymmetry-making magnetic field - as follows!^ 

^total ,asymm (0, (/), t) = Bshapeifj)) siri^O - Bstressit) COS^9 - Basymm{t) COS 0, (6) 

where Bshape{4>) and Bstressit) are given by the Equations ([2]) and ([3]), respectively, Basymmit) = MH(t), and 
H(t) is the asymmetry making magnetic field at time t. Note that it is possible to harness the asymmtery making 
field by intrinsic interface coupling between polarization and magnetization in multiferroic heterostructures.^^ 

When a sufficiently high stress (of right polarity making the product Xga negative) is generated on the 
magnetostrictive nanomagnet and the in-plane shape-anisotropic energy barrier is overcome, magnetization will 
rotate toward the in-plane hard axis (0 = 90°, (j) = ±90°). This is depicted in Fig. |BJ:. Magnetization initially 
can situate at any of the two degenerate stable states (0 = 0°, state 0 and 0 = 180°, state 1 with probabilities, 
Po = Pi = 0.5). The initial entropy of the system at time ti is —kJ^nPn ^n-Pn = kln{2). However, when stress 
is applied between times H and ^ 2 , the barrier separating the two stable states gets removed and the potential 
landscape becomes monostable in 0-space (at 0 = 90°) with entropy zero. A reduction of entropy k ln(2) must be 
dissipated as heat during this process, according to the Landauer’s principle. Note that the barrier is removed 
before applying any asymmetry-making field H. The rationale behind is to make the erasure process independent 
of barrier height, which determines the hold failure probability and also to resist the thermal fluctuations by 
making the monostable well deep enough. 

Between times t 2 and ta, an asymmetry-making field H is applied to deterministically rotate the magnetization 
to state 0 = 0°. This is depicted in Fig. [b}:. The asymmetry-making field tilts the potential landscape and the 
degree of this tilt should be sufficient enough to dissuade thermal fluctuations. Thereafter, both the stress and 
H are removed to complete the switching process. Note that the entropy of the system at time ^2 and onwards 
are zero. Hence there is no bound on minimum energy dissipation when magnetization traverses from time ^2 
to fa as the dissipation can be made arbitrarily small. Assuming (3/2)As(Tmaa; = —MH^ and cj) = ±90°, from 
Equation ([6]), the total energy becomes Etotai,asymm{9, t) = —{l/2)MHk sin^O — MH{t) cosO. We notice that as 
H goes from 0 to Hk , the minimum value of 0 goes from 0 = 90° to 0 = 0° continuously following the expression 
9min{t) = COS~^{H{t)/Hk). 

With adiabatic pulses between times ti and ta, the motion of magnetization is smooth and magnetization 
follows the minimum potential energy landscape. Therefore there is no lower bound of energy dissipation in the 
absence of thermal fluctuations. However, the Landauer’s principle remains intact since this is T = 0 K case 
and thus the energy dissipation proportional to temperature is also zero. If we incorporate random thermal 
fluctuations at finite temperatures, magnetization will get deflected uphill in the potential landscape even with 
pulses of very slow ramp and will incur energy dissipation, which is subjected to Landauer bound. 

It needs mention here that particularly in the presence of thermal fluctuations, magnetization may temporarily 
traverse on higher potential not only in-plane of the nanomagnet but also out-of-plane (i.e., when cj) ^ ±90°) and 
dissipate energy when it comes back to the lower potential. The simulation of magnetization dynamics results 
that a less than one degree of deflection in the out-of-plane direction can have an immense consequence. The 
key reason behind is that the out-of-plane demagnetization field is about a couple of orders of magnitude 
higher than Hk- When the magnetization gets deflected out-of-plane due to torque exerted on it, an additional 
torque of comparatively very high magnitude due to Hd comes into play, which makes the dynamics fast. It is 
true that this out-of-plane excursion causes power dissipation but switching also becomes fast, so that the net 
energy dissipation may be smaller compared to the case when it is assumed that magnetization is confined to 
the magnet’s plane. Since the Landauer limit of energy dissipation is very small, we should particularly take into 
account this significant effect due to the out-of-plane excursion of magnetization, considering realistic parameters 
for magnetization, e.g., magnetization damping a. 

The torque Te acting on the magnetization is derived from the gradient of potential landscape. Additionally, 
there is a random thermal field to incorporate thermal fluctuations.^^ We solve the stochastic Landau-Lifshitz- 
Gilbert equatiorP^^ of magnetization dynamics [Equation (jd])] and calculate the energy dissipation during the 
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Figure 7. Simulation results for Landauer limit of energy dissipation in a magnetostrictive nanomagnet. 

Distribution of energy dissipation at room-temperature (300 K) during the time interval t 2 — ti. It is important to note 
that for a few cases, the Landauer bound {kT ln{2) = 0.6932 kT) is violated but the mean of the energy dissipation 
(0.74 kT) exceeds the Landauer bound safeguarding the Landauer’s principle and the second law of thermodynamics. A 
moderately large (1000) number of simulations in the presence thermal fluctuations have been performed to generate this 
distribution. 

erasure of a bit of information. The energy dissipated in the nanomagnet due to Gilbert damping can be expressed 
as = fj Pd{t)dt, where r is the time taken during the erasure cycle [i.e., — ti in Fig. [BJo], and Pd{t) is the 

power dissipated at time t per unit volume given bj^ 

Pdit) = (7) 

Thermal field with mean zero does not cause any net energy dissipation but it causes variability in the energy 
dissipation by scuttling the trajectory of magnetization. 

2.2 Results 

The magnetostrictive nanomagnet is considered to be made of polycrystalline Galfenol (FeGa), which has the 
following material properties - magnetostrictive coefficient ((3/2)As): -1-150x10“®, saturation magnetization 
{Ms): 8x10® A/m, Young’s modulus (Y): 140 GPa, and Gilbert damping parameter {a): 0.025.1^1^ The 
dimensions of the nanomagnet is 100 nm x 90 nm x 6 nm, which ensures that the nanomagnet has a single 
ferromagnetic domain.I^I^ With the chosen dimensions, the Stoner-Wohlfarth switching field Hk becomes 
~0.01Ms. The values of stress (amax), strain (e), and asymmetric field {H) are 60.6 MPa, 433x10“®, and 0.01 
T, respectively. 

Figure [7] shows the distribution of energy dissipation at 300 K during the time interval t 2 — ti = 100 ns, 
when stress is ramped up from zero to the maximum value, which makes the potential landscape monostable 
and erases the bit of stored information. This decreases the entropy of the system by kln{2) = 0.6932 fc. 
A concomitant amount of energy must be dissipated according to Landauer’s principle. Stochastic violation 
of Landauer’s bound 0.6932 kT due to thermal fluctuations is observed but the mean energy dissipation does 
respect the Landauer’s bound. This signifies the generalized Landauer principle for small systems, which are 
prone to thermal fluctuations. The mean energy dissipation during the time interval — ^2 (when magnetization 
traverses from 9 = 90° towards 6 = 0°) is 0.07 kT, which does not have any bound since there is no entropy loss 
in the system. 













3. CONCLUSIONS 


Multiferroic straintronic devices turn out to be promising for energy-efficient computing in beyond Moore’s law 
era. Since stress anisotropy in the magnetostrictive nanomagnets contribute a symmetric term, it is necessary 
to have an asymmetric component to maintain the direction of switching. Intrinsic magnetization dynamics can 
contribute to such asymmetric term. With realistic parameters, the slight out-of-plane excursion of magnetization 
plays a crucial role in shaping the magnetization dynamics and therefore it is of paramount importance to 
consider complete three-dimensional potential landscape and solve full three-dimensional magnetization dynamics 
rather than assuming magnetization is always confined to magnet’s plane. Also intrinsic coupling between the 
polarization and magnetization can harness the asymmetric term. These straintronic devices can be utilized 
not only for digital computing but also for analog signal processing. We have shown as well that the ultimate 
limit (so-called Landauer limit) of energy dissipation is achievable in a magnetostrictive nanomagnet, linking 
information and thermodynamics. These hndings would hopefully stimulate experimental efforts to demonstrate 
that the magnetostrictive nanomagnets are suitable for exploring the thermodynamic limit of energy dissipation. 
The miniscule energy dissipation in these straintronic devices can be the basis of ultra-low-energy computing for 
our future information processing systems. This can open up as well some unprecedented applications that need 
to work with the energy harvested from the environment e.g., monitoring an epileptic patient’s brain to notify 
an impending seizure by drawing energy solely from the patient’s body. 
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